Evidence for Strain-Induced Local Conductance Modulations in Single-Layer 

Graphene on Si02 
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Graphene has emerged as an electronic material that is promising for device applications and for 
studying two-dimensional electron gases with relativistic dispersion near two Dirac points. Nonethe- 
less, deviations from Dirac-like spectroscopy have been widely reported with varying interpretations. 
Here we show evidence for strain-induced spatial modulations in the local conductance of single- 
layer graphene on Si02 substrates from scanning tunneling microscopic (STM) studies. We find 
that strained graphene exhibits parabolic, U-shaped conductance vs. bias voltage spectra rather 
than the V-shaped spectra expected for Dirac fermions, whereas V-shaped spectra are recovered 
in regions of relaxed graphene. Strain maps derived from the STM studies further reveal direct 
correlation with the local tunneling conductance. These results are attributed to a strain-induced 
frequency increase in the out-of-plane phonon mode that mediates the low-energy inelastic charge 
tunneling into graphene. 
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Since its experimental isolation in 2005,— graphene 
has emerged as an electronic material with superior 



propertiei 



1.2.3 



that are promising for device apphca- 



tionsi^'^'^ Additionally, graphene represents a unique 
opportunity to study two-dimensional electron gases^i^ 
with relativistic dispersion near two Dirac points in the 
Brillouin zone.- That is, the charge carriers in ideal 
graphene are massless Dirac fermions with a conical 
energy-momentum relationship, E = zszhkyp, where vp 
is the Fermi velocity, H = /i/(27r), h is the Planck 
constant, and k the momentumii Thus, the low en- 
ergy charge excitations of ideal graphene are expected 
to obey the Dirac equation rather than the Schrodinger 
wave equation^i In addition to the prediction of novel 
low-energy excitations, various interesting phenomena 
have been observed in bulk measurements of graphene, 
including the ambipolar electric field effectfi the inte- 
ger quantum Hall effect (IQHE) ^i^i^iiiS the presence of 
a minimum conductance 4e^/(fi.7r) in the limit of zero 
charge-carrierSj^i^ ultra high mobilities fi^ii^ and the ab- 
sence of quantum-interference magnetoresistance.-2iii Mi- 
croscopically, while a number of investigations have been 
conducted on single point spectroscopy of graphene us- 
ing scanning tunneling spectroscopy (STS),^ii^ii^iii in- 
formation on spatially resolved local density of states 
(LDOS) over extensive areas has been limited except for 
recent studies of epitaxial graphene on SiC(OOOl) and 
Ru(0001)M 

In contrast to the lack of spatially resolved spec- 
troscopic information, substantial high-resolution topo- 
graphic studies of graphene have been carried out using 
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STM. For instance, it is found that strong surface corru- 
gations, up to 1 nm variation in height on a lateral scale 
of 10 nm, appear in the topography of both graphene on 
Si02 substrate o^°i^^ and suspended graphene. These 
findings suggest that not only does the graphene devi- 
ate from a true two-dimensional system with a slight 
three-dimensional component, but the planar structure 
also deviates from a perfect honeycomb lattice. Such 
massive changes in the topology are expected to af- 
fect the LDOSi^i^^i^iiS^iS^ as well as the phonon modes 
of graphene. Recent findings of out-of-plane phonon- 
mediated inelastic tunneling from STS stuides of me- 
chanically exfoliated graphcne^^ have revealed the impor- 
tant role of phonons in the low-energy excitation spectra 
of graphene. Additionally, substrate-induced sublat- 
tice symmetry-breaking in the case of expitaxially grown 
graphene can give rise to energy gaps in graphene elec- 
tronic structures.^ Thus, surface corrugations and struc- 
tures associated with the substrate may induce signif- 
icant variations in the low-energy tunneling spectra of 
graphene. 

In this work we address the issue of possible effects of 
lattice distortions on the local tunneling conductance by 
STM studies of single-layer graphene on Si02 . Graphene 
samples were manufactured via mechanical exfoliation. 
An optical microscope was used to identify the location 
of single-layer graphene, and photolithographic processes 
were employed to attach gold electrical contacts to the 
graphene. Prior to STM measurements, the sample was 
annealed in a high oxygen environment at temperatures 
of 400 °C for 15 minutes to remove photoresist that re- 
mained from the photolithography processes; otherwise, 
atomic resolution would not have been achievable and 
spectroscopic measurements would have been contami- 
nated. The topographic and spectroscopic measurements 
were taken with a home built cryogenic STM, which was 
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capable of variable temperature control from room tem- 
perature to 6 K and was also compatible with magnetic 
fields. For studies reported in the following, we kept 
the measurement conditions at 77 K under high vacuum 
(< 10"'' torr) and in zero magnetic field. Topographic 
and spectroscopic measurements were performed simul- 
taneously at every location in a (128 x 128) pixel grid. 
At each pixel location, the tunnel junction was indepen- 
dently established so that the junction resistance of 4017 
was maintained across the sample. The differential con- 
ductance, [dl/dV), was calculated from the best polyno- 
mial fit of each current (/) vs. bias voltage {V) curve. 

Topographic STM measurements were made over nu- 
merous regions across the sample, as exemplified in 
Fig. 1(a) for a (2.2 nm x 5.0 nm) area. We found sur- 
face corrugations of up to ±0.5 nm over a lateral dis- 
tance of 10 nm, as shown by the histogram in Fig. 1(b), 
which is in good agreement with the results reported 
previously .^^'^^1^^ Additional confirmation of the sur- 
face corrugation was obtained from atomic force mi- 
croscopy (AFM) measurements over the same region, 
which recorded height variations of ±1.0 nm. Further to- 
pographic STM studies over a larger area of the graphene 
sample revealed similar landscape of surface corruga- 
tions, as shown in Fig. 1(c). The corresponding his- 
togram of the height variations is depicted in Fig. 1(d). 

To understand the physical origin for the graphene sur- 
face corrugations, we peformed AFM measurements on 
a bare Si02 substrate, which yielded height variations 
of ±2.0 nm. In Fig. 1(d) we compare the histogram of 
the graphene surface corrugations with that of the Si02 
substrate over a (8.6 nm x 8.6 nm) area. The apparent 
correlations between the two histograms suggest that the 
surface corrugations of the graphene sample result from 
the underlying roughness of the substrate. We further in- 
vestigated the possibility of graphene surface roughness 
resulting from the gold contacts on the graphene sam- 
ple, and found no apparent correlation between the sur- 
face roughness and the distance from the gold contacts. 
Additionally, extensive topographic surveys revealed no 
discernible atomic defects in our graphene sample. 

To quantify the lattice distortions of graphene on Si02, 
we performed fast Fourier transformation (FFT) of the 
topographical scan over the (2.2 nm x 5.0 nm) region 
in Fig. 1(a). The FFT in Fig. 2(a) revealed a distorted 
hexagon, indicating significant deformation in the lattice 
structure. Thus, spatially varying strain maps may be 
derived from studying the displacement fields u of local 
regions. Specifically, we may choose one of the principle 
lattice vectors as the a;-axis, and define the local dis- 
placement field \i{x, y) = UxX + Uyij as the difference 
of the local lattice vector from the equilibrium lattice 
vector. Following similar approaches previously devel- 
oped for electron microscopic studies of other materi- 
als,— we obtained the strain maps shown in Figs. 3(a)- 
3(d) for the sample area given in Fig. la: Fig. 3(a) il- 
lustrates the "scalar strain" map So{x,y) = |u(a;, y)/ao|, 
and Figs. 3(b)-3(d) correspond to the maps of the strain 




FIG. 1: Surface corrugations of graphene on Si02 substrate: 
(a) Topographic image of graphene surface over a (2.2 nm x 
5.0 nm) area, (b) Histogram of the graphene height variations 
over the area shown in Fig. 1(a). (c) Topographic image of 
graphene surface over a larger (8.6 nm x 8.6 nm) area, (d) 
Histogram of the graphene height variations over the larger 
area shown in Fig. 1(c) is compared with the histogram of 
a bare Si02 substrate over the same area, showing appar- 
ent correlation between the height variations of graphene and 
those of the underlying Si02 substrate, (e) Two-dimensional 
map of the Dirac voltage over the area in Fig. 1(a), showing 
relatively homogeneous Dirac voltage distributions, (f) His- 
togram of the Dirac voltage over the same area in Fig. 1(e). 



tensor components Sxx{x,y) = {dux/dx), Sxy{x,y) — 
Syx{x,y) = [(dux/dy) + {duy/dx)]/2, and Syy{x,y) = 
{duy/dy). From the strain maps we find that on the 
whole the upper section of the graphene in view was more 
strained while the lower section was more relaxed. 

In addition to the togographic images, spatially re- 
solved spectroscopic studies were carried out across the 
graphene sample to investigate possible correlation be- 
tween the surface corrugations and modulations of the 
LDOS. The differential conductance (dl/dV) as a func- 
tion of bias voltage (V) was obtained at each pixel of the 
sample area under investigation. We first noticed that at 
each pixel the corresponding spectrum always exhibited 
a minimum conductance at a finite bias voltage Vd ■ This 
finite Vd value, known as the Dirac voltage, corresponds 
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FIG. 2: Correlation of spatially varying topography and con- 
ductance maps of graphene: (a) FFT of the topographic map 
in Fig. 1(a), showing a distorted hexagon, (b) A conduc- 
tance ratio map taken at Vb = 0, revealing significant spatial 
variations in the conductance that are correlated with the to- 
pography in Fig. 1(a). Here the conductance ratio is defined 
as the local conductance relative to the mean conductance of 
the entire area in the field of view, (c) A conductance ratio 
map taken at Vb = 80 meV, showing weaker spatial variations 
in the conductance relative to Fig. 2(b). (d) A conductance 
ratio map taken at Vb ~ 240 meV, showing reduced spatial 
variations in the conductance relative to Figs. 2(b) and 2(c). 
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FIG. 3: Strain maps of graphene over a (2.2 nm x 5.0 nm) 
area: (a) A "scalar strain" map So{x, y), showing the spatially 
varying fraction of lattice distortion relative to the equilib- 
rium lattice constant ao = 2.46A. The three white dashed 
lines (upper, lower and vertical) indicate three different line- 
cuts of varying strain tensor components. (See Fig. 4.) (b) 
A map for the strain tensor component Sxx{x,y). (c) A map 
for the strain tensor component Sxy{x,y). (d) A map for the 
strain tensor component Syy{x,y). 



to an offset in Fermi energy, and the predominant mecha- 
nism for this offset may be attributed to charge impurities 
on the Si02 surface. For the area shown in Fig. 1(a), the 
corresponding two-dimensional map and the histogram 
of the Dirac voltage Vd(j:, j/) are illustrated in Figs. 1(e) 
and 1(f), respectively, showing relatively uniform vari- 
ations without apparent reference to the topography in 
Fig. 1(a). In contrast, the spatial map of tunneling con- 
ductance at Vb = {V — Vd) = 0, Fig. 2(b), reveals strong 
correlation with the topography in Fig. 1(a). As Vb in- 
creases, the spatial variations in conductance across the 
sample decrease steadily, as exemplified by Fig. 2(c) for 
Vb = 80 meV and Fig. 2(d) for Vb = 240 meV. The lack 
of correlation between the Vd{x, y) map in Fig. 1(e) and 
the low-energy conductance map in Fig. 2(b) suggests 
that charge impurities are unlikely the primary cause for 
the spatial modulations in the LDOS. 

To understand possible correlations between the spa- 
tial modulations of tunneling conductance and the local 
strain, we further examined the tunneling spectra and the 
strain tensor components along three different line-cuts 



(upper, vertical, lower) across the sample, with the three 
line-cuts indicated by the white dashed lines in Fig. 3(a). 
The spatial evolution of the spectra associated with the 
upper hue-cut is given in Fig. 4(a), which are mostly 
parabolic (U-shaped) around Vb = 0, as exemplified by 
the solid curve in Fig. 4(b). We further note that the 
spatial distribution of these U-shaped spectra correlates 
well with the region of strained graphene according to 
the maps in Figs. 3(b)-3(d). 

If we employ the scenario of out-of-planc phonon- 
mediated inelastic tunneling^l by allowing the out-of- 
plane phonon frequency Wa as a fitting parameter and 
taking T = 77 K, as elaborated in the next paragraph, 
we find that the theoretical fitting (dashed curve in 
Fig. 4(b)) agrees reasonably well with the experimental 
spectra for a phonon frequency huja = 44 meV. In con- 
trast, the spectra of the lower line-cut in Fig. 4(c) are 
more Dirac-like (V-shaped) around Vb — 0, as exempli- 
fied by the solid curve in Fig. 4(d). Similarly, the region 
where such V-shaped spectra dominate coincides with 
the area of relaxed graphene with nearly vanishing strain 
tensor components as shown in Figs. 3(b)-3(d). More- 
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FIG. 4: Correlation of the tunneling spectral characteristics 
at T = 77 K with the strain tensor components: (a) A series of 
differential conductance {dl /dV) vs. bias voltage {Vb) spec- 
tra along the upper line-cut in Fig. 3(a), showing primarily 
parabolic U-shaped spectra around Vb = 0. (b) A repre- 
sentative U-shaped normalized conductance {dl /dV)/{I /V) 
spectrum (solid curve) in the strained region. A theoreti- 
cal curve with an out-of-plane phonon frequency fkJa = 44 
meV — is shown by the dashed line. The inset shows the 
theoretical curves at T = 0. (c) A series of {dl /dV)-vs.- 
Vb spectra along the bottom line-cut in Fig. 1(a), showing 
primarily Dirac-like V-shaped spectra around Vb = 0. (d) A 
representative V-shaped normalized conductance spectrum in 
the relaxed region together with the theoretical fitting curve 
(dashed line) with a phonon frequency huja = 26 meV. The 
theoretical curve at T = is illustrated in the inset, (e) A 
series of {dl / dV)-YS.-VB spectra along the vertical line-cut in 
Fig. 3(a), showing evolution from primarily U-shaped spectra 
to V-shaped spectra, (f) Comparison of two representative 
normalized conductance spectra from strained to relaxed re- 
gions. The theoretical fitting curves correspond to a higher 
phonon frequency hcua = 41 meV associated with the strained 
region and a smaller frequency htOa = 24 meV with the re- 
laxed region. 

over, the theoretical fitting to the spectrum in Fig. 4(d) 
yields a smaller out-of-plane phonon frequency hua — 26 
meV relative to that in the strained region. For the ver- 
tical line-cut that moves from a region of strained to re- 
laxed graphene, the corresponding spectra in Fig. 4(e) 
exhibits evolution from spectra with larger huja values 
to those with smaller fkUa values while the strain tensor 



components decrease according to Figs. 3(b)-3(d). Direct 
comparison of two representative spectra from strained 
to relaxed regions with a decreasing out-of-plane phonon 
frequency is shown in Fig. 4(f). The apparent correlation 
of spectral characteristics with the strain maps together 
with their lack of correlation with the Diract voltage 
Voix^y) suggest that strain-induced conductance mod- 
ulations and phonon-mediated inelastic tunneling rather 
than charge impurities are the primary cause for spectral 
deviations from Dirac behavior in graphene. 

Next, we speciiy our fitting to the tunneling conduc- 
tance in Fig. 4 based on the three-band model that con- 
siders the mixing of nearly-free electron bands at the zone 
center F with the Dirac fermions at the zone edge K and 
K' through coupling with the out-of-plane phonons.^^ Ac- 
cording to this scenario, the tunneling DOS from STM 
measurements is dominated by the DOS of the nearly- 
free electron bands, iVr(w). In the limit of |w| ^ W 
where W « 6eV denotes the Dirac electron bandwidth, 
Nr{ijj) in the absence of any gate voltage (/i = 0) is given 

S'l' 1 {uj + iS) 
Tr\uj — ha — 2ji^i(a; + «o)| 

Here « 3.3eV, and the real and imaginary parts of 
the electron self-energy Ei^i(a; -I- i5) = S'i_i(w -I- i5) + 
iT,'l i{uu + iS) satisfy the following energy dependence)^ 

S'i^i(c. + iS) cx (c. - COa) log l^^l -I- (c. -f c.,) log 

T,i i{oj + iS) CC ~Q{\uj\ ~ LUa)\uJ — Sgn(^Uj)LL)a\- (2) 

In our data analysis we assume Ua as a fitting parameter, 
and also allow a small constant zero-bias offset in Eq. (1), 
which may be attributed to enhanced LDOS due to im- 
puritiesii^i^ For T = 77 K, we replace uj by iuJn and sum 
over the fermion Matsubara frequencies a;„ with a pref- 
actor {1/(3) where /9~^ = 77 K. Given that the thermal 
energy at 77 K is much smaller than hiOa, however, the 
thermal smearing effect is insignificant, as manifested by 
Figs. 4(b) and 4(d). 

Finally, we note that the hiUa values derived from the 
phonon-mediated tunneling scenario appear to vary sig- 
nificantly not only within different regions of a substrate 
but even more among different substrates, ranging from 
24 ^ 44 meV found in our experiments to 60 meV and 
100 meV reported in RefsJ^ andi^, respectively. In the 
case of graphene grown on SiC substratesfi^ superlattice 
modulations are known to have induced modifications to 
the electronic bandstructures,^^ which may account in 
part for the significantly larger gap features. Although 
we cannot pinpoint the exact origin of gap variations for 
gaphene prepared under different conditions, the signifi- 
cant strain-induced effects appear to be at least partially 
responsible for the differences. On the other hand, it 
seems unlikely for the out-of-plane phonon mode alone 
to account for the significant variations of the empirical 
huJa values among the tunneling spectra of graphene on 
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different substrates. Further investigation appears nec- 
essary to address this issue. 

In summary, our findings of substrate-induced lattice 
strains in graphene and the apparent correlation between 
the strain tensor components with modulations in the 
local tunneling conductance provide new insights into 
mechanisms that can strongly influence the low-energy 
excitations of Dirac fermions. From the viewpoint of de- 
vice applications of graphene, our studies suggest that 
strain-induced variations in the electronic properties are 
limited to relative low excitation energies, around a few 
tens of meV's. Therefore, device applications involving 
finite gate or bias voltages on the order of eV's2.!ii£ will 
not be sensitive to the substrate-induced strain effect. On 
the other hand, it is conceivable to apply strain-induced 



effects to controlling the low-bias electronic properties of 
graphene-based devices. In this context, our findings of 
strain-induced energy gap variations in graphene can be 
of technological relevance. 
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